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Abstract 

Electrons on a square lattice with half a flux quantum per plaquette are 
considered. An effective description for the current loops is given by a two- 
dimensional Dirac theory with random mass. It is shown that the conductivity 
and the localization length can be calculated from a product of Dirac Green's 
functions with the same frequency. This implies that the derealization of 
electrons in a magnetic field is due to a critical point in a phase with a 
spontaneously broken discrete symmetry. The estimation of the localization 
length is performed for a generalized model with N fermion levels using a 
1/iV-expansion and the Schwarz inequality. An argument for the existence of 

two Hall transition points is given in terms of percolation theory. 
PACS numbers: 71.55J, 73.20D, 73.20J 
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Lattice models for the two-dimensional electron gas in a strong magnetic field are of 
increasing interest because of recent experiments with GaAs heterostructures Jl| . Moreover, 
lattice models offer new ways to the study of the quantum Hall transition. The purpose of 
this article is to discuss a new concept for the investigation of the localization length for the 
lattice model. 

The tight-binding Hamiltonian for non-interacting electrons on a square lattice with 
magnetic flux reads in Landau gauge 

H = -t ^[e 2i ^/^ c +(r + e x )c(r) + c + (r + e„)c(r) + h.c] + £ V(r)c + (r)c(r). (1) 

r r 

e x>y are lattice unit vectors, and c + and c are fermion creation and annihilation operators, 

respectively. V(r) is a potential which represents an additional structure or disorder on 

the lattice. The dispersion for V = and half a flux quantum per plaquette (0 = <f>o/2) 

is E = ±2^1/ cos 2 ki + cos 2 k 2 . A linear (relativistic) approximation for E ~ around the 

four nodes kj = ±7r/2 is possible. For H this is still very complicated ||, and a further 

simplification is useful to lift the degeneracy of the four nodes. This can be achieved by the 

introduction of a next nearest neighbor hopping term and a staggered potential in H ||,|J. 

This modification is motivated by the network model |J. Expansion of k = (ir/2, vr/2) + p 

for small p vectors leads to the large scale approximation by the Dirac Hamiltonian H' = 

a ■ p + (T3M with Pauli matrices <jj. Disorder, originally introduced in H by a staggered 

random potential V, is described in H' by a random mass M |§. Physical quantities can 

be obtained from the Green's function 

G 12 \ (iuo + M h Y 1 

\G 2 i G22) \ h + iu-M) 
where the Fourier component of h is p\ — ip 2 - For instance, the frequency- dependent con- 
ductivity is given by Kubo's formula 

e 2 

(Txx(w) = — ^ 2 ^r 2 (G ji /(r,0;ia;)G' J v i (0, r; -i^)), (3) 

where Gjj'(r, 0; iu>)Gjij(0, r; — iuS) = \Gjf(r, 0; iu)\ 2 . The localization length ^ is defined as 
the decay of the function Cjji(r,u) = (\Gjj>(r,0;iu)\ 2 ) in space. There exists a relation of 
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the above expression, which is composed of Green's functions at frequencies with opposite 
sign (retarded and advanced Green's functions), with a product of Green's functions at the 
same frequency. This follows from the block matrices in (2) 

Gjj(r,r'; -iu) = -Gj> jr (r',r;iu), G jf (r,r';-iuj) = -Gjj>(r',r;iw) (j ^ j). (4) 

This identity reflects the Lorentz-covariance of the Dirac theory It implies \Gjj(r, r'; ioj)\ 2 = 
—Gjj(r,r';iuj)Gj'jr(r,r';iuj) and \Gjf(r, r'\ iu))\ 2 = —Gj>j(r,r';iu)Gjji(r,r';iu). That means 
only the Green's functions with one frequency is required for the evaluation of transport or 
localization properties in the present model. 

The physics of the electrons can be understood as the statistics of current loops, created 
by the magnetic field. Depending on the potential there are local current loops with two dif- 
ferent directions. This observation is central for the understanding of the Hall transition dis- 
cussed subsequently. The current loops are composed from the Dirac fermions by the creation 
and annihilation of particle-hole pairs. The direction can be reversed globally by a time- 
reversal transformation. In terms of the Green's function this means G(u) — > —a^G{uj)a^. 
In general the system is not invariant under this transformation because one direction of the 
currents is favored. The favored direction characterizes the electronic state of a Hall plateau. 
The transition between the Hall plateaux corresponds to a symmetric point with uj = M = 
which is also a critical point with a divergent localization length. The discrete symmetry 
at this critical point plays a fundamental role for the lattice electrons in a magnetic field. 
Taking now a random mass, which is symmetrically distributed around zero, this symmetry 
is spontaneously broken for the average Green's function ||. It will be discussed in the 
following that the localization length also diverges at two critical points in the presence of 
disorder where the symmetry is broken. This is an extension of a previous work where the 
Hall conductivity was calculated. Critical points with a divergent localization length in a 
symmetry broken regime are known from Anderson localization in systems without magnetic 
field and dimensionality d > 2. However, the physics of the latter is different because there 
are no current loops. As a manifestation of this difference a continuous symmetry appears || 
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instead of discrete symmetry found for localization in the presence of a magnetic field. The 
large scale approximation for the continuous symmetry leads to a non-linear sigma model 
contrary to the Dirac theory in the presence of the magnetic field. 
The Schwarz inequality can be applied to get a lower bound for C: 

\(G jf ,{r,0;iu)G fjl »(0,r;iu))\ 2 < (\G jf ,(r, 0; iu)\ 2 )(\G jlf »{0, r; iuj)\ 2 ) (5) 

with j" = j, f = f or f = f, f" = j. Writing C'^r, u) = (G^(r, 0; ico)G ff (0, r, iu)) and 
C'jj,(r,uj) = (Gjf(r,0;iuj)Gjij(0,r,iuj)) one has |C'y(r,a;)| < \Cjf(r,uj)\. C will be used 
subsequently because it easier to calculate than C. The average correlation functions are 
translational invariant. Therefore, the corresponding Fourier components C(k,u) can be 
used to calculate the localization length 
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(6) 

k=uj=0 y ' 



A simple calculation for the pure system gives = |M| _1 . To evaluate the correlation 
function C'(r,u) a generalization of the Hamiltonian H' is introduced which describes iV 
levels of fermions per site 0: H aa ' = H$ a ' - 5M? a '(r 3 (a, a' = 1,2,...,N) with Hff* = 
(a ■ p + {M)as)5 aa ' . The distribution of the real symmetric random matrix (a Gaussian 
Orthogonal Ensemble |]) 5M is given by (5M? a ' 5Mf a "') = (g/N)5 aa "' 5 a ' a " 5 r y. That 
means only random fluctuations couple the N different Dirac systems. The product of two 
Green's functions with the same frequency can be expressed as a functional integral [ 10 1 



Gjj„ (r, r'; w)Gf$ (r", r'"; ico) = - f ^v,^^,,,^^ exp(-S'i) II d$ r d$ r (7) 

J r 

with the action (sum convention for the level index a) 

Si = z<7 w [-($, (H + icuao)®) + E SMr'i^r • ^)\- (8) 

r 

cr w = sign(a;) and the field is = (^"j, The first component is Grassmann and the 
second is complex. The complex component is added to normalize the functional integral 
in (7). It also provides a transparent representation of the product of Green's functions 



with the same frequency as required for C and C . Averaging with Gaussian distributed 
fluctuations yields (exp(— Si)) = exp(— S 2 ) with 

S 2 = -ia u (a>, (H + iu* )$) + -| ■ <7 3 i>?')« ■ a 3 ^). (9) 

This action can be mapped onto another effective field theory [[ll|]: The interaction in S 2 
can also be created by complex 2 x 2-matrix fields Q, P and a complex Grassmann field 
ij) which couple to the composite field J2a=i &r^r- The new field does not depend on a, 
the level index of the N level fermions. The level degree of freedom can be eliminated by 
integrating out the field $ in the functional integral. This leads to exp(— NS(Q, P, ip)) with 



the action [[11 



S(Q, P, V) = - £[Tr 2 (Q r0 r 3 ) 2 + Tr 2 (P r a 3 ) 2 ] + log det(# - 2Q + iua ) 

9 r 

1 4 r 

-logdet(F + 2iP + zW ) + - ^ J {\)^{k)^ k ^. ky d 2 k + 0{N- 2 ). (10) 

/ t ,/i'=i 

Since the number of fermion levels A^ appears in front of the new action, the limit A^ — > oo 
corresponds to a saddle point integration for the fields Q and P. I(k) is the matrix of the 
Gaussian fluctuations around the saddle point solution Q = Q + 5Q and P = iQ + 5P, 
where Qo is the A^ — ► oo solution Q = —(l/2)[iri<T + M s <r 3 }. M s is a shift of the average 
Dirac mass and rj shifts the frequency in the Green's function. Introducing m = (M) + M s 



as the effective (renormalized) Dirac mass the imaginary shift is rj = o^\J e~ 2lT / 9 — m 2 for 
\m\ < m c = e~ n / g and zero otherwise ||11|| . 

The details of the derivation of S(Q,P,ip) from S 2 and the evaluation of saddle point 
integration can be found in Ref. |]TT| . However, (10) can also be reconstructed using the fact 
that / exp(— Si)d... = 1. The A^ — > cxo-limit gives Qo = —iPo which implies that the large 
A^-terms cancel each other in the action (10). First order fluctuations do not contribute 
at the saddle point. And second order fluctuations appear as quadratic forms of 5Q and 
5P with the matrix I. Therefore, the condition Jexp(— NS(Q, P, ip))d... = 1 is satisfied 
if a quadratic form of I appears with a complex Grassmann field. This result reflects the 
super symmetric construction of the functional integral. 



I(k = 0) 



For the asymptotic behavior of a large correlation length it is sufficient to consider k ~ 0. 

If k = the matrix I reads 

fl/g-2afx 2 2/5 \ 

2/g + 4a\n\ 2 

2/g + Aa\fi\ 2 

V 2/3 l/g-2afi* 2 J 

where // = m + irj, a = J{\fi\ 2 + k 2 )- 2 d 2 k/47c 2 and (3 = J k 2 (\fi\ 2 + k 2 )- 2 d 2 k/An 2 . The 
inner 2x2 matrix of I(k = 0) is diagonal with positive matrix elements 2/g + 4a|/i| 2 . It 
leads always to a finite correlation length. However, this part of I does not contribute to C 



because Cjy 



T 1 



TTj] . Therefore, only the projected submatrix (P14 is the 
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(12) 



projector on the indices 1,4) 

l/ 5 -2/i 2 a + ] u 2 4fc 2 + 

2/3 + Bxk 2 l/g-2/j* 2 a + y? 2 A x k 2 

with positive constants A\ y Bi is of interest here. It describes the critical behavior of the 

model due to one vanishing eigenvalue for k = at two critical points m = ±m c . The 

critical points are also characterized by a vanishing imaginary part of the saddle point Qo- 

The Fourier components of C are of the form C'(k,u = 0) = {2(3 + bk 2 )/(a' + b'k 2 ), 

where b and b' are finite constants and a' = \l/g — 2afi 2 \ 2 — 4/5 2 . Equ. (6) yields for the 

correlation length £ of C 



i oc 2{b'/a' - 6/2/3) 1/2 ~ 2^yjd. 



(13) 



Because b' is always finite the divergence of the correlation length is determined by 1/af. 
The behavior of a' distinguishes the regime of \m\ < e~ n ^ 9 and \m\ > e~ n ^ 9 . This leads to 



f ~ £o(l/0+ -log|m|)- 1/2 for \m\ > e^ /g 

7T 



(14a) 



and 



i ~ £o(e~ 2 " /9 - m 2 )" 1/2 for |m| < e^ /s . (146) 

Thus £(m) diverges linearly at m = ±e _7r / 9 (then 77 = according to the saddle point 
solution), and it is finite elsewhere. The existence of two critical points, instead of the single 



6 



critical point in the perturbative approach of Ref. is a consequence of the spontaneous 
symmetry breaking due to disorder which creates a non-vanishing density of states between 
the particle and hole band of the Dirac theory. This non-perturbative phenomenon was 
discussed in Ref. @,|],[n] . 



The correlation function C cannot be calculated directly from S(Q,P,ip) because it not 
expressable as a simple form in the large N limit. Therefore, only C can be used as a lower 
bound for the decay of C due to inequality (5). As a consequence, £ bounds the localization 
length £j from below. An alternative approach, based on a different effective field theory 
with N — 1, indicates also a square root divergent localization length [12|. The exponent 



v fa 1 for the localization length was found experimentally for the metal-insulator transitions 



(i.e. near the mobility edge) in AlGaAs/GaAs heterostructures Q and in Si MOSFETs [13 

On the other hand, numerical results for network models and lowest Landau level ap- 
proaches f5,14-|T6f suggest that the localization length £j diverges at the critical point with 



an exponent v « 7/3. 

The existence of two transition points is in disagreement with the single transition point of 
other approaches to the integer Hall transition |5|JT^-|i7| . This is either due to different mod- 
els (Landau level approach versus lattice fermions) or due to the fact that the two transition 
points were not seen because of limitations of the numerical resolution. (The distance of the 
two transitions is 2e~ n ^ 9 which is extremely small for weak disorder. It can be hidden by 
finite size effects.) Two transition points are also supported by the percolation picture for 
the lattice model. This can be seen considering a lattice where the space-dependent Dirac 
mass M r is either zero with probability p or M with probability 1 —p. The correlation length 
£ is infinite for p — 1 because M = is the symmetric point. For p < 1 not all lattice sites 
have a zero mass. However, there is typically an infinite cluster of lattice sites with zero mass 
for p c < p < 1 according to percolation theory ||18|| . The correlation length of the classical 
percolation system is finite for p < p c and diverges like £ ~ £o(Pc — vY with the critical 
exponent v = 4/3. The mean value of M r at the critical point p c is (M) = (1 —p c )M = M c . 
It seems plausible that its value is non-zero even if a renormalization ((M) — > m) by quan- 



turn fluctuations of the fermions is taken into account. Application of the transformation 
G — > a 3 Ga 3 to the Green's function G(u = 0) replaces M by —M. Thus — M c is also a 
critical point. 

In conclusion, it was shown that the transport of electrons on a square lattice in a strong 
magnetic field and random potential can be described by a product of Dirac Green's function 
with the same frequency. Therefore, the physics is characterized by a discrete symmetry. 
The localization length diverges at two different critical points. The critical exponent was 
estimated by v > 1/2. 
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